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abstract. The present study aimed to determine the complete citrate synthase (git A) and heat-shock protein (groEL) gene sequences of 
Anaplasma bovis and to infer phylogenetic relationships within the genus Anaplasma. Multiple alignments from single and concatenated 
sequences of the 16S rRNA, git A and groEL genes of the genus Anaplasma were subjected to phylogenetic analyses. Percent identities of 
A. bovis nucleotide sequences were found highest with A. phagocytophilum in gltA (65.4%) and groEL (79.8%). Single gene phylogenetic 
tree results assumed similar phylogenetic positions within the genus Anaplasma, except for A. bovis. However, consensus and concatenated 
sequence phylogenetic trees showed similar results, revealing 2 subgroups within the genus. 
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The genus Anaplasma currently recognizes 6 distinct spe- 
cies: A. phagocytophilum, A. platys, A. marginale, A. cen- 
trale, A. ovis and A. bovis. The reclassification was mainly 
based on the phylogenetic information derived from 16S 
rRNA (complete representation of all member species) and 
heat-shock operon or groEL (only selected member species) 
gene data sequences [6]. A recent study identified a poten- 
tially novel Anaplasma sp. in Japan (herein provisionally 
referred to as Anaplasma sp. Japan), which revealed phylo- 
genetic divergence in the 16S rRNA, gltA and groEL genes 
from any recognized Anaplasma spp. [21]. However, the 
previous studies present paucity of information on whether 
the use of secondary structures was employed in their phy- 
logenetic analyses, which is known to produce better tree 
resolution [15]. In addition, the widely accepted 16S rRNA 
gene based phylogenies are sometimes inconsistent [16], 
which is probably due to the propensity of the 16S rRNA 
gene to recombination/horizontal or lateral gene transfer 
phenomenon [1]. Therefore, other genes like the gltA [12] 
and groEL [11, 13] can be alternatively used to clarify phy- 
logenetic relationships. 

A. bovis infects circulating monocytes [5]. This particular 
species has been mainly analyzed only using the 16S rRNA 
gene [6]. Previous phylogenetic analyses of the genus Ana- 
plasma using the groEL gene did not include yet A. bovis 
[6, 21] due to the unavailability of the data sequence during 
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the time of analyses. The present study generally aimed to 
molecularly characterize and analyze A. bovis based on gltA 
and groEL genes and to infer phylogenetic relationships 
within the genus Anaplasma using individual and multi-lo- 
cus approach (including the 16S rRNA gene). Phylogenetic 
analyses were performed with or without the consideration 
of secondary structures, using maximum likelihood (ML) 
and Bayesian Inference (BI) methods. 

Blood sample from a feral raccoon (Procyon lotor) [18] in 
Hokkaido, herein referred to as R499, was used. The sample 
was previously tested to be 16S rRNA-positive for A. bovis 
(1,387 bp; GenBank accession number GU937020) and was 
stored at -30°C. The DNA was extracted and stored as previ- 
ously described [21]. 

The designing of primers, determination of the partial gltA 
and groEL sequences of A. bovis by PCR, genome walking 
and DNA sequencing strategies were performed as described 
previously [21]. Primers used in the present study are shown 
in Table 1 . The negative control used was double distilled 
water. Instead of using an A. bovis DNA, the positive control 
used was A. platys. 

The gltA and groEL sequences were translated into de- 
duced amino acids (dAA) and were manually trimmed to 
include only the sequence of interest (generally from the 
start to stop codon). Percent identities were computed as 
previously described [21]. Multiple sequence alignments 
(MSA) were performed as suggested by Hall [7] or by using 
PROMALS3D [15], which considers secondary structures 
for protein coding genes. Subsequent analyses with and 
without using the secondary structure information were 
performed using raxmlGUI [19] by general time reversible 
(GTR) model. Analyses by ML with prior best model testing 
using MEGA 5 [20] and by BI using MrBayes 3.2 [17] were 
also employed. For the protein coding genes, ML analyses 
were performed using MEGA 5 with prior best model test- 
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Table 1 . Oligonucleotide sequences of primers used in this study 



Primer name Oligonucleotide (5'— >3') Reference 

PCR amplification for a partial gltA sequence 
CS7F2 ATGR*TAGAAAAW*GCTGTTTT [21] 

HG1085R ACTATACCK*GAGTAAAAGTC [10] 

Fib GATCATGAR*CAR*AATGCTTC [9] 

AnaCS1085Rl ACTATACCK*GAGTAAAAR*TC This study 

PCR amplification for a partial graEZsequence 
EEgrolF GAGTTCGACGGTAAGAAGTTCA [3] 

Anagro712R CCGCGATCAAACTGCATACC [21] 

Anagrol22F AAATACGGTW*GTCACGGG This study 

Anagro649R CTTTCTTCR*ACAGTTATAAG This study 

Genome walker gene-specific primers (gltA) 

ABgl-46Rl AATGCAGCTGCTCCCGCACTTAAGCAAGTA This study 

ABgl-lR2 GTTGAGCCCACCATTCTTACTGTAGATGTA This study 

ABgl-307Fl CAATATAGCGATCGCTATGGAAGAAATAGC This study 

ABgl-338 F2 TTGCAGGATGAATACTTCATAGAGAGGAAG This study 

Genome walker gene-specific primers (groEL) 

ABgr-361Fl TAAAGGCAAAAGAGGCTGTTCTTACGG This study 

ABgr-385F2 CGGC TCTT ATGT CCAT GAGA CGTG AA This study 

ABgr-1010GRF3 ACAGTGCATCTTCCAGCATAGAAAGTAG This study 

ABgr- 1 1 1 7F4 AGCTTTCAGGTGGTGTGGCTGTGCTGAAAG This study 

ABgr-1132F5 TGGCTGTGCTGAAAGTTGGTGGATCAAGTGA This study 

ABgr-96Rl TTAGGTCCAGCAGTACATCCAACAGCA This study 

ABgr-30R2 GGATGTGTACTATCTCCCTAATGGCTTA This study 

ABgr-25R3 TCTCCCTAATGGCTTTATCCAACTTCTC This study 

*Degenerate primers: R=A or G, W=A or T, K=G or T 



ing, while BI was performed using MrBayes 3.2 guided by 
the prior best model test results from MEGA 5. All analy- 
ses utilizing MEGA 5 were estimated using 100 bootstrap 
replications. Selected representative sequences from species 
which had available information on the 3 different genes 
were concatenated (herein referred to as "supermatrix"). Su- 
permatrices from sequences were analyzed using MrBayes 
3.2. Tree results from MrBayes and raxmlGUI were viewed 
using the FigTree vl.3.1 (http://tree.bio.ed.ac.uk/software/ 
figtree/). Phylogenetic tree inputs were generated from each 
MSA (with secondary structure guidance) that were concat- 
enated for the supermatrix by BI and were used to construct 
of a majority-rule consensus tree using the Dendrosope [9]. 
Using the SplitsTree4 [8], the supermatrix was subjected to 
phylogenetic network analysis (by NeighborNet method) 
and was tested for the presence of recombination phenom- 
enon using the PHI test. 

Initially, partial gltA (382 bp) and groEL (385 bp) sequenc- 
es of A. bovis from R499 were obtained. After the genome 
walking procedure and sequencing, the complete gltA (1,239 
bp; JN588561) and groEL (1,644 bp; JN588562) sequences 
were determined. In comparison with other Anaplasma spe- 
cies, the gltA and groEL sequences of A. bovis were found 
closest to A. phagocytophilum (Table 2). Within the genus 
Anaplasma, species shared at least 60.8 and 60.4% identi- 
ties in the nucleotide and dAA sequences of the gltA gene, 
respectively and 77.2 and 89.4% identities in the nucleotide 
and dAA of the groEL gene, respectively (data not shown). 



In the 16S rRNA phylogenetic analyses, 2 subclades were 
seen: (1) a subclade containing A marginale, A. centrale and 
A. ovis and (2) a subclade containing A. phagocytophilum, A. 
platys, A. bovis and Anaplasma sp. Japan (Fig. 1). A. bovis 
also frequently formed a cluster with Anaplasma sp. Japan. 
In the gltA phylogenetic analyses (Fig. 2), topologies re- 
vealed the 2 subclades observed in the 16S rRNA trees. In 
the groEL phylogenetic analyses, some positions within the 
genus Anaplasma changed depending on whether nucleotide 
or dAA sequences were used, but the 2 subclades were still 
frequently observed (Fig. 3). On the other hand, trees gener- 
ated from the supermatrix also revealed the 2 subclades (Fig. 
4). 

A. bovis consistently formed a cluster with A. phagocyto- 
philum, A. platys and Anaplasma sp. Japan in the 16S rRNA 
phylogenetic trees. This finding varied from the tree results 
of Dumler et al. [6], in which their 16S rRNA phylogenetic 
analysis placed A bovis closer to A. centrale and A. ovis, but 
was similar to that of Ooshiro et al. [14] and Doan et al. [4], 
in which A. bovis formed a cluster with A. phagocytophilum 
and A. platys. 

For the gltA and groEL gene phylogenetic analyses, the 
subclade groupings of the different taxa appear to be consis- 
tent when protein secondary structures were considered in 
the MSA construction, than when nucleotide sequences were 
used. The graEL-based trees generated in the present study 
also varied from those of Dumler et al. [6] as sequences 
of A. bovis, A. ovis, A. centrale and A. platys were not yet 
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Table 2. Percent identities of Anaplasma bovis from R499 with other Anaplasma spp. based on 16S rRNA, 
gltA and groEL gene sequences 









% Identity* 




Organism 


16S rRNA 




gltA 


groEL 




Nucleotide 


Nucleotide Amino Acids 


Nucleotide Amino Acids 


Anaplasma phagocytophilum 


96.7 


65.4 


62 .4 


79.8 91.4 








NC007797 




Anaplasma platys 


97 


60.8 


61.1 


79.6 92.9 




AY077619 


AY077620 


AF478129 


Anaplasma ovis 


95.5 






77.4 90.9 




AF309865 






AF441131 


Anaplasma marginale 


95.4 


74.6 


62.4 


77.5 90.9 








NC012026 




Anaplasma centrale 


95.5 


73.6 


60.4 


77.2 90.9 








NC013532 




Anaplasma sp. Japan 


96.5 


63.7 


61.6 


79 89.4 




JN055357 




JN55361 


JN55359 



*Gaps were not considered in the computation. Note: the accession numbers of the sequences compared are 
indicated below the cells. 



79 



100 



23 



99 



' Rickettsia prowazekii AJ235272 



Anaplasma marginale str. Florida NC012026 
65 1 Anaplasma ovis AF309865 
100 \r-Anaplasma ovis AF318945 

LAnaplasma centrale NC013532 
63 100 Anaplasma bovis R499 GU937020 
L Anaplasma bovis U03775 



33 



64 



Anaplasma sp. AY570538 
[ Anaplasma sp. AY5705389 
Anaplasma sp. Japan AB196721 
Anaplasma sp. Japan EU368727 



69L Anaplasma sp. Japan JN055357 

[Anaplasma phagocytophiium "equi" M73223 
[ 46 

Anaplasma phagocytophilum str.HZ NC007797 

[■ Anaplasma platys AY077619 

901 

1 Anaplasma sp. BomPastor AF318023 

'Ehrlichia ruminantium NC006831 
Ehrlichia muris GU358691 



Ehrlichia chaffeensis NC007799 
34 l|10Q Ehrlichia canis CP000107 

7fH Ehrlichia ovina AF318946 
Ehrlichia ewingii U96436 

Wolbachia pipientis AE017196 



Neorickettsia risticii NC013009 



0.04 

Fig. 1. Phylogenetic trees based on the 16S rRNA gene. Analyses were performed by maximum likelihood 
(general time reversible model) using 100 bootstrap replications employed in raxmlGUI [18]. Rickettsia 
prowazekii was set as the outgroup. 



included in their analyses. The groEL sequences of A. cen- 
trale, A. ovis [13] and A. platys [12] were only determined 
at a later time. Dumler et al. [6] pointed out the ambiguities 
among Anaplasma spp. and the arbitrary position of A. bovis 
within the Anaplasma species clade in the various phyloge- 



netic analyses they performed. 

Comparing the single gene or the multi-loci phylogenetic 
trees, the consistently observed result was the formation of 
the 2 subclades when secondary structures were considered. 
Moreover, the resulting topologies corroborated with our 
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Anaplasma marginale Florida NC012026 
Anaplasmo marginale Maries NC004842 
Anaplasma centrale Israel NC013532 
'Anaplasma centrale Aomori AF304141 
68 1 Anaplasma phagocytophilum HZ NC007797 
Anaplasma phagocytophilum Webster AF304136 
Anaplasma phagocytophilum Equi AF304137 
i Anaplasma sp. Japan JN055362 
' Anaplasma sp. Japan JN055361 
| Anaplasma platys AY530807 



100 



100 'Anaplasma platys DQ525686 



■ Anaplasma bovis R499 JN588561 



Ehrlichia ruminantium NC006831 
Ehrlichia ewingii DQ365879 
Ehrlichia can/s_CP000107 
Ehrlichia muris AF304144 
100 1 Ehrlichia chaffeensis NC007799 



I Ehrlichic 
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■ Wolbachia pipientis AE017196 



• Neorickettsia risticii NC013009 



— Rickettsia prowazekii AJ235273 



0.2 



Fig. 2. Phylogenetic trees based on gltA with consideration of the protein secondary structures. Analyses were 
performed by the Bayesian method (Jones-Taylor-Thornton model) employed in MrBayes 3.2 [16]. Values 
in the nodes represent posterior probability values expressed in percent. Rickettsia prowazekii was set as the 
outgroup. 
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Anaplasma marginale Florida NC012026 
Anaplasma marginale Maries NC004842 
Anaplasma centrale Israel NC013532 
Anaplasma centrale AF414867 
Anaplasma ovis AF441131 

Anaplasma ovis FJ460441 
Anaplasma phagocytophilum HZ NC007797 
Anaplasma sp. Japan from Ixodes persulcatus JQ622144 
Anaplasma phagocytophilum Horse AF172162 

Anaplasma phagocytophilum HGE Agent AF172163 
92 

Anaplasma sp. Japan JN055360 
96 Anaplasma sp. Japan JN055359 

Anaplasma sp. Haemaphysalis douglasii JQ685509 
Anaplasma platys AF478129 
L- Anaplasma bovis R499 JN588562 
— Ehrlichia ruminantium DQ647012 

Ehrlichia ewingii AF195273 
■ Ehrlichia cants U96731 
Ehrlichia sp.AB454077 
Ehrlichia chaffeensis L10917 
Ehrlichia muris AF210459 

Wolbachia pipientis AE017196 



-Neorickettsia risticii NC013009 



Rickettsia prowazekii AJ235272 



0.4 



Fig. 3. Phylogenetic trees based on groEL genes with consideration of the protein secondary structures. Analyses 
were performed by the Bayesian method (Jones-Taylor-Thornton model) employed in MrBayes 3.2 [16]. Values in 
the nodes represent posterior probability values expressed in percent. Rickettsia prowazekii was set as the outgroup. 
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Fig. 4. Results of the concatenation approach by the Bayesian method using MrBayes 3.2 [16]. All nodes in the figure 
revealed posterior probability values of 1 or 100%. Rickettsia prowazekii was set as the outgroup. 



previous findings [20], in which the Anaplasma sp. Japan 
was found to be a potentially novel species. The absence of 
statistical evidence of recombination (using PHI test) and the 
subsequent result of the phylogenetic network analysis (by 
NeighborNet method) on the concatenated alignment also 
supported the reliability of the tree results. PHI tests are used 
to test MSAs for the presence of recombination, which can 
obscure the results of phylogenetic analyses [2]. 

The present study documented the first molecular analy- 
ses of A. bovis based on complete groEL and gltA gene 
sequences and inferred phylogenetic relationships within the 
genus Anaplasma with the inclusion of new sequence data. 
Results clarified the phylogenetic position of A. bovis and 
established the existence of 2 subclades within the genus 
Anaplasma. This information can serve as a guide to future 
phylogenetic studies using the same genus. 
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